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Abstract. In this survey, our aim is to emphasize the main known limitations 
to the use of Wigner measures for Schrodinger equations. After a short review 
of successful applications of Wigner measures to study the semi-classical limit 
of solutions to Schrodinger equations, we list some examples where Wigner 
measures cannot be a good tool to describe high frequency limits. Typically, 
the Wigner measures may not capture effects which are not negligible at the 
pointwise level, or the propagation of Wigner measures may be an ill-posed 
problem. In the latter situation, two families of functions may have the same 
Wigner measures at some initial time, but different Wigner measures for a 
larger time. In the case of systems, this difficulty can partially be avoided by 
considering more refined Wigner measures such as two-scale Wigner measures; 
however, we give examples of situations where this quadratic approach fails. 
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1. Introduction 

In this survey, we briefly review some successful applications of the Wigner 
measures for classical limits of Schrodingcr equations, and discuss in more detail 
the limitations of this tool. Although most of the material presented is essentially 
known, we feel that a survey presenting the dis-advantages of Wigner measures in 
a clear unified picture is timely and useful. 

Wigner measures are a very valuable tool for describing high frequency and ho- 
mogenization limits for oscillatory PDEs, possibly with periodic coefficients. The 
Wigner measure is a phase space measure that allows to describe weak limits of 
quadratic quantities (the observables) of a (solution) family of functions which only 
converges weakly itself. The basic idea goes back to E. Wigner who used such a 
phase space approach in quantum mechanics for semi-classical approximations in 
1932 [68]. In the 90's, Wigner functions and their limiting measures aroused the 
interest of mathematicians in the USA (e.g. [63]) and, independently, in Europe. 
As a variant of L 2 defect measures (see [66, 38]) such objects where used as a tech- 
nicality in proofs in the frame of the analysis of crgodic properties of eigenfunctions 
for the Dirichlet problem (see [64, 69, 22, 49]), with more systematic studies of 
semiclassical measures by P. Gerard and E. Leichtnam (see [37, 41], or the survey 
of N. Burq, [4]). The term "Wigner measure" was used first in the French work 
"Sur les mesures de Wigner" of P.-L. Lions and T. Paul [54]. 

Adaptations to the case of Schrodingcr operators with periodic coefficients 
and applications of the method to general problems were given by P. Gerard, 
P. Markowich, N.J. Mauser and the late F. Poupaud [37], [58], notably their joint 
paper [42] , where a general theory of the use of Wigner measures for homogeniza- 
tion limits of energy densities for several wide classes of dispersive linear PDEs is 
laid out. In the context of our work we also want to mention the more recent use 
of Wigner measure for proving resolvent estimates ([5, 50, 51, 35]) following an 
idea of proof by contradiction of [53]. 

The method of Wigner measures allows to treat some (weakly) nonlinear equa- 
tions, for example the semiclassical limit of the coupled Schrodinger-Poisson sys- 
tem, first done in '93 in [54] and [57], both works using "smoothed Wigner func- 
tions" as a technical step to prove the non-negativity of the limiting measure and 
both crucially depending on the use of mixed states — an assumption that could 
only be lifted in ID so far [73]. The case of the inclusion of the additional difficulty 
of a periodic crystal potential was solved in [3] , where the general theory of Wigner 
series, as Wigner measures in the context of a Bloch decomposition of L 2 , was laid 
out. 

In general however, Wigner measure methods are not suitable for treating non- 
linear problems and even their use for linear problems has severe limitations. The 
main aim of this paper is to list some (rather explicit) examples where the use 
of Wigner measures is not appropriate and thus unveil in a clear and concise way 
the inherent strengths and shortcomings of the Wigner measure approach for the 
semiclassical limit of time-dependent Schrodinger equations. 



1.1. Setting of the problem. The semi-classical limit of the Schrodinger 
equation can be seen as a model problem for the kind of homogenization limits 
studied by this method. The rigorous mathematical development of Wigner mea- 
sures was motivated by this problem. By semi-classical limit we mean the limit of 
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the (scaled) Planck constant tending to in the Schrodingcr equation, which reads 

(1.1) ied t ip £ = -—^ e + V(x)i> e x£R d ,t£R 

(1.2) Vf t =o = V>f, t€R. 

Here e stands for the scaled Planck constant, tp e = tp e (t, x) £ C is the wave function, 
and V(x) £ R is a given potential. The wave function may be vector-valued and the 
potential then is an Hermitian matrix: in this situation, (1.1) is a system. Clearly, 
this limit is a high frequency limit and can only exist in some weak sense. 
From the viewpoint of physics, the main interest is not on i\) e itself, but in quantities 
which arc quadratic expressions in ip e , e. g. the position density 

(1.3) n s (t,x) = \r(t,x)\ 2 . 

Of course, this quadratic operation docs not commute with the weak e —> limit 
of t/j s . The importance of getting limits of quadratic quantities like (1.3) and the 
prevalence of the distinct scale e of oscillations in the solutions of (1.1) make Wigncr 
measures the "good" tool for this problem and generally for the homogenization of 
energy densities of time dependent PDEs with a scale of oscillations. 
The Wigner transform of a function / £ L 2 (R d ) is defined by 

(1.4) w*[f](x,0 = (2ir)- d J j{ x - £ l)y( x + e l) C " 1<d7] - 

In this definition, e is an arbitrarily introduced parameter, the scale of the Wigncr 
transform. w £ [f](x,!;) is real-valued, but in general not positive. Now let w £ (t,x^) 
be the Wigner transform of ip e (t, x), then for n e (t, x) defined by (1.3) it holds that 

(1.5) n E (t,x)= w E (t,x,£)d4, x£R d 

Jmf 

Other quadratic quantities in ip E can be obtained by taking (higher) moments in 
the ^-variable of w £ (t, x, £). In the case of vector- valued /, the Wigner transform 
is matrix- valued: the product of / and / is replaced by a tensor product in (1.4). 
One still has (1.5) provided one takes the trace of w £ . 

In order to study the convergence properties of w E , in [54] the following space 
was introduced: 

A = & £ C (Ri x Rf) \(F c <p)(x,ri) G L 1 (R d - C (R d x )) } , 
with 11^^11^(0=/ supl^Kx,?/)^, 

Jr™ x 

where the Fourier transform is defined as follows: 

(1.6) = = 

It is easy to verify that A is an algebra of test functions and a separable Banach al- 
gebra containing S{R d x R d ). This space immediately allows for a uniform estimate 
on the Wigner function : 

Proposition 1.1. Let -0 e be a sequence uniformly bounded in L 2 (R d ). Then 
the sequence of Wigner transforms w E [tjj E ] is uniformly bounded in A' . 

It follows that, after selection of a subsequence, 

(1.7) w E m(x,0 -±w°(x,t) in A'. 

It can be shown that w (a;,£) is a non-negative measure on the phase space: 
the semi-classical or Wigner measure of the sequence ip e , which is not necessar- 
ily unique. Note that the Wigner transform w e is in general real, but may also 
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have negative values, whereas the limit w° is non-negative, thus justifying the term 
"Wigner measure" . In the vector- valued situation, one uses matrices of test func- 
tions and w° is a non- negative hermitian matrix of measures, which means that 
w = (idJ •) with non-negative measure and absolutely continuous with 
respect to w% and • . 

We could work also with uniform bounds and convergence in the distribution 
space S as the dual of the Schwartz space, as it was done e.g. in [42] without 
recourse to the space A . 

In the scalar case, when applying the Wigner transform to equation (1.1), a 
kinetic transport equation will result, the so-called "Wigner equation". The formal 
limit under e — > of this transport equation leads to the following Vlasov equation 
for the Wigner measure w°: 

(1.8) d t w° + £ ■ V x w° - V x V(x) ■ V 6 w° = 

If this limit can be made rigorous, and the above equation is well posed as an initial 
value problem, the Wigner measure oiip e (t, x) at a positive time can be obtained by 
solving the equation with the initial Wigner measure as data: the Wigner measure 
is constant along the Hamiltonian trajectories (a;(i, y, 77), £(i, y, ??)) 

j d t x(t,y,rj) = £(t,y t r)) ; x(0,y,T))=y, 

\d t Z{t,y,r)) = -V x V(x(t,y,r))) ; £(0, y, r,) = V . 

For describing the limit in the position density n £ , the following definition will 
be necessary: 

Definition 1.2. A sequence {f £ } 6 uniformly bounded in L 2 , is called e- oscillatory 
if, for every continuous and compactly supported function (j> on M. d , 



(1.9) lim / 0/^(0 



d£ — ► 0. 

R — >oo 



This definition is presented here in the form that is used in [42], following 

^ 2 

[37] and [41]; in [54] the equivalent condition that l/e d f £ (^/s) is a relatively 

compact sequence in Ai(M. d ) is used instead. Heuristically this means that the 
wavelength of oscillations of f e is at least e. A sufficient condition for (1.9) is 

3k > such that e K D K f £ is uniformly bounded in L 2 oc . 

If ij} £ is e-oscillatory, it is possible to pass to the limit in (1.5), and we find 

(1.10) / w°(x,d£) = (w-) limn e (a;) = n°(x) in V, 

If now the Wigner measure at a positive time t > can be obtained in an appropri- 
ate sense by the evolution equation (1.8), this identity serves to get the desired limit 
in the density. In the context of quantum mechanics the strength of the method 
lies in the fact that the transport equation of the Wigner measure (1.8) is identical 
to the transport equation from classical statistical physics. So according to the 
"correspondence principle" the quantum problem converges to a classical problem 
in the limit where the "quantum parameter" e vanishes. 

The following diagram gives a schematic sketch of the method used to obtain 
homogenization limits of quadratic quantities by Wigner measures. 
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Wigncr transform 



(weak limit) 



.Transport equation 



^Time evolution 
(*D0) 



r(t,x) 



n s (t,x) = \il> e {t,x){< 



e 

(weak limit) 



w°(t,x,d£) 



?? 



n°(t,x) 



If it is possible to do the operations on the left hand side of this diagram, the 
limit of n £ can be found by applying (1.5) to the Wigncr measure at a positive time 
t > 0. So it is necessary 

• to have a unique Wigner measure of the data, 

• that the transport equation for w° is a well-posed problem such that a 
solution exists up to a relevant time t > 0, 

• the ^-integral (zero-th moment) of w° at t > must be equal to the limit 
n°(t,x). 

For the third step to be valid, tp e (t, x) has to be e-oscillatory as stated above, 
which means that there must be no oscillations at faster scales than e. If the e- 
oscillatory property is imposed on the data, it will be preserved by (1.1) for positive 
times. The second step, existence and uniqueness of a global solution to the trans- 
port equation of the Wigner measure, is known to hold for a large class of linear 
scalar problems, but it turns out to be much more complicated for systems (see 
Section 3); however in nonlinear settings this point is a mostly open question. 

Note that in the nonlinear case the non-uniqueness of the Wigner measure of 
the sequence of solutions ip E (t, x) somewhat corresponds to the non- uniqueness of 
the (weak) solution w°(t,x) of the limiting nonlinear Vlasov equation; clearly, we 
are not able to pick one particular solution of the ill-posed PDE by a semiclassical 
limit of the unique solution of the nonlinear Schrodinger equation. 



1.2. Some successful applications of Wigner measures: scalar case. 

For e-independent potentials, as in (1.1), the convergence of the Wigner function to 
a Wigner measure as a solution of a Vlasov equation was systematically described 
in [54]. If the potential V is sufficiently smooth (C 1,1 ), then the Wigner measure 
for (V ,6 ) £ is given in terms of the Hamiltonian flow associated to |£| 2 /2 + V(x). 
On the other hand, if V is not sufficiently smooth, uniqueness for the Hamiltonian 
flow fails, and the above mentioned convergence is not guaranteed; see [54] for an 
exhaustive list of examples and results. 

Note that actually the first mathematical study of "semiclassical measures" 
for linear Schrodinger equations was done by P. Gerard [37] in the context of the 
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setting in a crystal, where Bloch waves are considered. 

2 

(1.11) ied t r = - £ -M° + V r (J) ^ ; ^(0, x) = Vo ( J) , 

where Vr is lattice-periodic. In the independent later work [58] the names "Wigncr 
Bloch functions" and "Wigner series" for their limits were coined, phase space 
objects that are essentially obtained via the definition (1.4) by replacing the Fourier 
integral by a Fourier sum, keeping the position variable x in whole space and 
restricting the kinetic variable to the torus (the Brillouin zone) as the dual of the 
lattice. 

Wigner measures have also proven successful in the study limits of quadratic 
quantities of the Schrodingcr-Poisson system: 

(1.12) ied t r =-yAf + W ; VfU = V'f, 

for a potential V^ L = V^^'ip 6 ] 2 ) which depends on n £ = \ip £ \ 2 . The main problem 
in applying the Wigner measure method to this equation is the low regularity of 
w° , which generally is only in A! . This means that the transport equation for w° 
can only be fulfilled in a rather weak sense, so in general a nonlinear expression in 
w° will have no meaning. 

The applications of Wigner measures to nonlinear problems which are available 
so far hold for Hartree type nonlinearities, where the nonlinear potential V^ L is 
given by 

(1.13) Vf, L = J U(x-y)n £ (y)dy 

for some suitable U and n e given by (1.3). For the Schrodingcr-Poisson system, 
one usually considers U(x) = l/|x| when d = 3 (and in the absence of background 
ions). The semiclassical limit results for nonlinear cases in [54] and [57] are pos- 
sible only for the case of a so-called "mixed state", i.e. by considering infinitely 
many Schrodinger equations. In this case convergence can take place in a stronger 
sense since a uniform L 2 -bound on the initial Wigner transform can be imposed by 
choosing very particular initial states (see [59] for a comprehensive discussion). 

Note however that even in the case of such strong assumptions that is possible 
only if we replace the one pure state Schrodinger equation by infinitely many mixed 
state Schrodinger equations, the results are not strong enough for obtaining a limit 
that satisfies the conditions for ensuring unique classical solutions of the Vlasov- 
Poisson system (see [75]). The theory of global weak/strong solutions of the Vlasov 
equations is laid out, for example, in [26, 55] and [18, 36]), with famous non- 
uniqueness results notably for measure valued initial data in [74]. 

For the Schrodingcr-Poisson system 

.2 



" 2 



ied t ^ e = -— Arl> e + V£ L ip e 



with b = b(x) G L 1 (IR d ), the Wigner measure associated to ip e was computed in 
[72] for WKB-type initial data 



ij £ (0,x) = ^1 pl{x)e lMx)/e . 



It was proven that the Wigner measure is given in terms of the solution to an 
Euler-Poisson system, before the solution to the latter develops a singularity. We 
will see in §2 that in this special case, more can be said on the semi-classical limit 
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of tp £ . In particular, its pointwise behavior can be described, showing phenomena 
that the Wigner measures do not capture. 

In the special case of only one space dimension, a global in time description of 
the Wigner measure without any mixed state setting and the corresponding strong 
assumptions on the initial data is given [73]. Consider 



(1.15) 



£ 2 

led^ 6 = - — d^ E + V% L ilf, t>0, 
V E | t =o = V'f, x£ R, 



in the case where b £ ^ is the mollification of a function be L 1 P\L 2 , and is the 
mollification (with the same mollifier) of a sequence ip e bounded in L 2 (R). Then 
the Wigner transform of ip £ converges to a weak solution of the Vlasov-Poisson 
system 



(1.16) 



d t f + £d x f-Ed ( f = 0, f\ t =o = fi- 
d x E = b(x) - f m fd£, i£l, t > 0. 



Following [56] we denote by weak solution on the interval [0, T] any pair (E, /) 
consisting of E G (BV n L°°)([0, T] xl) and / G L°°(R+, M + (R 2 )) such that 

(1) \/<f> e C r °°([0,T] x R 2 ), 3 q4> g BV([0,T] x R), 

<t>{t, Of{t, x, d£ = d x g (j> . 

(2) E(t,x) = f* 00 (b(y)-f R f{t,y,t)d£)dy a. e, 

(3) V^€C c °°((0,T)xR 2 ), 

hf + <f>xtf)dx dt- f ( E ( <t> i f{dC)dxdt = d, 



E(t,x) 



10 JR 2 JO 

where E(t,x) is Vol'pert's symmetric average: 

E(t, x), if E is approximately continuous at (t, x), 

\{Ei(t, x) + E r (t, x)) if E has a jump at (t, x), 

where Ei(t,x), E r (t,x) denote the left and right limits of E at (t, x). 
(4) 3s > 0, / G C°^([0,T),^ c s (R 2 ) and f(0,x,0 = ^(x.f) in i^ c s (R 2 ). 
We refer to [73] for precise statement of the weak convergence result. We point out 
that the existence of weak solutions is proved in [7 4] , but the difficulty is that there 
is no uniqueness of these solutions: a counterexample is given in [56]. We close the 
section by pointing out that this latter result applies in particular to initial data 
which are e-independent or of the WKB form (f E (x) = p(x)e lS ( x ^ e , with p G L 2 (R). 

1.3. Some successful applications of Wigner measures: the case of 
systems. In the nonlinear case, several results are available in the case of the 
Schrodinger-Poisson system 



(1.17) 



where the total position density is defined as 

oo 

(i.i8) p £ = E A Mi 2 

j'=i 





j G N, x G R d , t > 0, 


V£ L = U*p s , 




i>j\t=o = 0j, 


j G N, x G R d , 
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with X € j E R_|_, Y^'jLi = 1- The Wigner transform for mixed states is defined, in 
analogy to (1.4), as 

(1.19) w*(t,x,0 = J^yJ Rd e - ^^ (t,x+^L,x f ) dy, 
where z £ (i, r, s) is the so-called mixed state density matrix defined by 

oo 

(1.20) z e (i,r, S ) = ^A^(t,r)V^(M), r,seR d . 

i=i 

Note that z £ (t, r, s) is the integral kernel of the density operator p e in L 2 , which is 
trace class with Ti(p E ) = YlJLi = 1- For more details we refer to [54, 57]. A 
crucial property of the mixed state is the fact that a uniform L 2 -bound on wj(x, £) 
holds if 

(1.21) ? a 

3=1 

This makes it possible to improve the sense of convergence for w £ (t, x, £) to weak-L 2 . 

In [3], the case of the Schrodinger-Poisson system was considered, with an 
extra Bloch potential (as in (1.11)). Note however that the analysis does not allow 
band crossings. Indeed, the analysis in terms of Wigner measures of systems face 
several difficulties in presence of eigenvalue crossings. 

The analysis of [42] for systems cover the case of matrix valued potentials with 
eigenvalues of constant multiplicity in (1.1). Let us denote by Xj the eigenvalues of 
V and by ILj the associated projectors, 1 ^ j ^ N. Then any Wigner measure w° of 

N 

ip e decomposes as w° = ^^u> 0j where the measures uP^ satisfy TLjiv^ILj — u> 0j 

and transport equations in the distribution sense 

d t w^ + £ ■ V x w°' j - V X A,- ■ V cW ° J - [Fj(x, 0,w°' j ] 

(the matrix Fj depends on Ilj, Xj, 1 ^ j ^ N, see [42, 43] for precise formula). As 
soon as the eigenvalues are not of constant multiplicities, this analysis is no longer 
valid and there may happen energy transfers between the modes which cannot be 
calculated in terms of Wigner measures. Such a phenomenon has been well-known 
since the works of Landau and Zener in the 30's (see [52] and [70]). It has been 
discussed from Wigner measures point of view in the articles [33] and [29] for a 
larger class of systems. Eigenvalue crossings appear for Schrodinger systems in 
the frame of quantum chemistry where one can find a large variety of potentials 
presenting such features (see [27]). In Section 3, for 

V(x)=( Xl X2 ), 
\ x 2 -xi J 1 

we explain how one can modify Wigner measure so that a quadratic approach still 
is possible and give examples where this approach fails. 

This article is organized as follows: we first describe in Section 2 the limitation 
of Wigner measures in scalar situations for e depending potential or nonlinear 
situations; then in Section 3, we focus on the case of systems. 

2. Limitations of Wigner measures in the scalar case 

The Wigner measure method is rather limited when nonlinear problems are 
treated, as can already be seen from the discussion of the above results. Note 
however that some limitations are present even in the linear case. We list below four 
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families of problems for which the loss of information due to the Wigner measure 
analysis is rather serious. 

2.1. Ill-posedness in the linear case: wave packets. We first discuss a 
situation where this approach fails in a linear setting, by recalling a case where the 
Cauchy problem (1.8) is ill-posed. Such an example for scalar Schrodinger equation 
is given by F. Nier [61]. In that case, the potential is e-dependent, as in (1.11), but 
is decaying instead of lattice-periodic. The problem is that the Wigner measure 
does not contain enough information on the properties of concentration of wave 
packets. Thus, Nier introduces in [62] a larger phase space, and a refined Wigner 
transform which takes into account the spread of the wave packets around the 
point of concentration involved. Similar ideas can be found in [60] and [28]. The 
quadratic approach does not fail once it is refined. The same kind of difficulties 
appear in systems with matrix-valued potentials presenting eigenvalues crossings 
(see Section 3). 

As a particular case of [61, 62], consider the problem: 



(2.1) = --Af + C/yf, 

(2-2) ^(0,,) = ^ (f), 

where U is a short range potential and uq S Ran W- = Ran W+, the wave operators 
to the classical Hamiltonian — ^A + U: uo = W-J-^ 1 ^'- = W+T^ip^. Introducing 
the solution u to 

id t u = -^Au + U (x) u ; U| t=0 = Uo, 

we see that: 

... . 1 ftx 
By assumption, we have u(t, x) ~ e l i A J 7 ~ 1 (ip±) as t — > ±oo, which implies 



2t as t — > ±oo, 



(it)"/ 2 

where we recall that the Fourier transform is normalized like in (1.6). Back to the 
initial unknown function ip e , this yields: 

Proposition 2.1 (from [61, 62]). Letip 6 be the solution to (2.1)-(2.2), where 
U is a short range potential and uq € RanW_ = Ran W + with uq = W-T~ x i\)- = 
W+ T~ Then the Wigner measure of the (pure) family ip £ is given by 

1 I*- (7 ) 



w°(t,x,() = \ "J" 



2 

cfo;<g> ift<0, 



2 



dx®St=* ift>0. 



|t|* 

As mentioned in [62], given two functions ip_ and ip'_ such that |^_| = \ip!_ |, 
one should not expect |^ + | = \^' + \, even in space dimension one. This is a first hint 
that the propagation of Wigner measures is an ill-posed problem in this context. 
This argument is made more precise in [61] (this example was not resumed in the 
complete paper [62]). In space dimension one, assume that the potential U is even 
U(x) = U(—x), and fix T > 0. Let ip £ be a solution to (2.1) such that its Wigner 
measure satisfies 

w°(-T,x,g) = 5 x =- Xo <g> <5 C=?0 , 

w°(T, x, = (1 - R 2 )Sx= X0 ® S i=(a + i? 2 4=-x ® . 
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where R = |-R(£o)| = \R{— £o)| is the reflection coefBcient related to the scattering 
operator S = W^_W-. Define V> E (i,x) = ip e {— 1,~ x )- Then since U is even, ip £ 
solves the same equation as ip e , and its Wigner measure satisfies 

w°(-T, x, = (1 - R 2 )S X =- X0 <8 S s=So + R 2 5 X=X0 ® <y c =_ Co , 
w; (T, x, £) = 5 x=3;o (8) <%=| . 

Define -0 e = V'f + V'l where the ijj^s solve (2.1) and whose initial Wigner measures 
are such that: 

w%(-T,x,0 = (l-R 2 )5 x= - X0 ®6^ ; w° 2 (-T, x, $ = R 2 S x=Xo <g> <5 5= _ ?0 . 

Then w° (-!>•,£) = 55° (-T, x, £), and unless = or 1, i<j (T,x,£) 7^ w°(T,x,£): 

^(T,x,£) = ((l- J R 2 ) 2 + J R 4 )5 :c = :Co ®^=£ +2ii 2 (l- J R 2 )5 :c= _ :Co ®^ = _ €o . 

Therefore, the Cauchy problem for the propagation of Wigner measures is ill-posed 
in this case: knowing the Wigner measure at time t = —T does not suffice to 
determine it at time t = +T. 

2.2. Caustic crossing and ill-posedness in a nonlinear case. We now 

give another example on a nonlinear problem, taken from [7], which may be viewed 
as a nonlinear counterpart of the above example. Consider a Schrodinger equation 
with power-like nonlinearity: 

(2.3) isd t r = ~&r+z 2 m i,d r ■, r\ t=0 = ao(x)e-^ , 

where x € M. d , d 1. The asymptotic behavior of tp e is given in [7] for any time. 
Note that other powers in the nonlinearity are also considered, provided that the 
power of £ in front of the nonlinearity is well chosen. 

PROPOSITION 2.2 ([7]). Let a e 5(R d ), and consider ip e solution to (2.3). 
Then the asymptotic behavior of tp E is given by \\ip 6 (t) — v £ (t)\\ L 2( R d- ) — > as e — ► 0, 
where: 



«o ( -A- ) e'^tt-D ift<l, 



v £ (t,x) 



(i-ty/z v \i-t 



Za [ — — ) e' 2 ^'- 1 ' ift > 1. 



(t - iyi* u \t-l, 

Here Z = T o S o T , where T stands for the Fourier transform (1.6), and S 
denotes the scattering operator associated to idtU-j- ^Au = |u| 4 / d u (see e.g. [19]). 

Like in Subsection 2.3, we infer the Wigner measure of tp E : 

2 

I x \ 



1 



\t-l\ d 
1 



Za 



( - \ 















dx <8> S^ = -2L- if t < 1, 
2 

dx (g> ^=-2- if t > 1. 



Unless I ao 1 2 = |-Zao| 2 , w° has a jump at the caustic crossing. The Wigner measure 
solves the transport equation with a singular source term: 

d t w°+^ V x w° =4=0® (\Za Q (C)\ 2 - |ao(OI 2 )#®<*t=i- 

The pathology is even more serious: the following result was established in [8] in 
the one-dimensional setting, and its proof extends to any space dimension. 
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Proposition 2.3. Let d > 1. 

(1) There exists ag G <S(R ) such that the Wigner measure w° associated to ip e 
solving (2.3) is discontinuous at t = 1: 

lim w°(t,dx,d£) ^ lim w (t,dx,d£_). 
t-»i- t-»i+ 

(2) There exist two (pure) families (V'|)o<e<i) j — 1,2, solutions to (2.3) (wiift 
different initial profiles aoj), whose Wigner measures are such that iv® = 
for t < 1 and iuj 7^ uj" /or t > 1 . 

Sketch of the proof. The point is to show that one can find (ao,j)j=i,2 
such that |ao,i| = |ao,2 1 and |Z(ao,i)| 7^ |-Z(ao,2)|- Since very few properties of the 
scattering operator S are available, the first two terms of its asymptotic expansion 
near the origin are computed, following the approach of [40] (the first term is 
naturally the identity). In the case of L 2 -critical nonlinear Schrodingcr equation 
considered here, we have, for ip- G L 2 (M. d ) and < 5 <C 1, 

= H- - ^ 1+4/d / u (-t) (|C/ (i)V-| 4/d £M*)V>- ) dt 

J — 00 

+G L2(Rd) (V+ 8 / d ) , 

where Uo(t) = e^ A . The proof of the above identity relies on Strichartz estimates 
and a bootstrap argument; a complete proof can be found in [16]. The idea is then 
to consider 

o ,i = a G S(R d ), a , 2 = a e ih , with h G C°° (K d ;R) . 
We proceed as in [8] . Denote 

/+00 
U (-t)(\U (t)iP-\ 4 / d U (t)^)dt. 
-OO 

Obviously, 

\To S W-)| 2 = S 2 U„| 2 + 2<5 2+4 / d Re (vD 5 ^) + (>+ 8 / d ) , 

and we have to prove that we can find ip_ G 5(M d ), and h G C 00 (M rf ; K), such that 

Rc(^T^( J PV-)) ^Rc(FWv)^(P(V^))) =--R{i>-,h), 
where ^ is defined by 

If this was not true, then for every ^_ G 5(R d ), the differential of the map n 1— ► 
R(if)—,h) would be zero at every smooth, real- valued function h. An elementary 
but tedious computation shows that 

D h R{i>_,0)(h)^0, 

with h(x) = \x\ 2 /2 and tp-(x) = e - ' 1 '"/ 2 . The computations uses the fact that the 
evolution of Gaussian functions under the action of the free Schrodinger group can 
be computed explicitly. We refer to [13] for more detailed computations. □ 

Remark 2.4. A similar result can be established in the case of the Hartree 
equation with harmonic potential studied in [15, Sect. 5]. For 7 > 1 and oq G 
l ), consider: 



x 



2 



(2.4) ied t ^ + -^ = l^^ + ei(\x\-T*m')^ ; VjU = «o 
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Like in [9] in the case of a power-like nonlinearity, the harmonic potential causes 
focusing at the origin periodically in time; outside the foci, the nonlinearity is negli- 
gible, and near the foci, the harmonic potential is negligible while the nonlinearity is 
not. Like for Proposition 2.2, its influence is described in average by the scattering 
operator associated to idtu + ^Au = 7 * |w| 2 ) u, whose existence was proven 
in [44, 48]. Following the same approach as in [8], one can prove the analogue of 
Proposition 2.3 in the case of Eq. (2.4). 



ied t i/j E = ~ £ —A^ e + V(t,x)ip e + eF E (t,x)ip e ; 4> e \ t=Q = a {x)e lMx)/E , 



2.3. WKB analysis and weak perturbations. Consider a Schrodinger 
equation with an 0(e) perturbation and a WKB data: 

f! 

2 

where V and F £ are real- valued. If we assume that F e has an expansion as e — > 
of the form F e = i*b + £-Fi + . . ., then a formal WKB analysis yields ip e « ae l ^l e , 
where: 

M+l\V<f>\ 2 + V(t,x) = ; 0| t=o = 0o- 

dta + V0 • Va + -aAcf) = —iFoa ; a | t=0 = a o- 

The first equation is the eikonal equation. Classically, it is solved locally in space 
and time, provided that V and 0o are smooth, by considering the Hamiltonian flow 

d t x(t,y) = £(t,y) ; x(0,y)=y, 

d t Z(t,y) = -V x V(t,x(t,y)) ; £(0, y) = Vfoiv)- 



(2.5) 



See e.g. [25, 46]. In addition, if cj>o and V are subquadratic, then one can find a 
local existence time which is uniform with respect to x G M. d : 

Assumption 2.5. We assume that the potential and the initial phase are smooth 
and subquadratic: 

• V G C°°(Rt x and d%V G L<g c (R t ; L°°(R d x )) as soon as \a\ > 2. 

• (f> a G C°°(R d ), and d£(f> Q G L°°(R d ) as soon as \a\ > 2. 

Lemma 2.6 (from [12]). Under Assumption 2.5, there exist T > and a unique 
solution cik G C*°°([0,T] x M d ) to: 

dk|t=0 



(2.6) d t eik + -|V0 ei k| 2 + ^(M) = 



This solution is subquadratic: <9"</> e ik G L°°([0,T] x R d ) as soon as \a\ ^ 2. 

Essentially, for ^ t ^ T, the map y i— > x(t, y) given by (2.5) is a diffcomor- 
phism of R d : for ^ t ^ T, the Jacobi determinant 

J t (y) = detVj,a;(t, y) 

is bounded away from zero. For ^ t T, no caustic is formed yet, and the second 
equation is a transport equation. It is an ordinary differential equation along the 
rays of geometric optics: introduce A given by 



A(t,y) :=a{t,x{t,y))y/Jt^). 

The transport equation is then equivalent, for ^ t ^ T, to: 

(2.7) d t A(t,y) = -iF Q (t,x(t,y))A(t,y); A(0, y) = a {y). 

We note that since F is real-valued, the modulus of A is independent of t G [0, T] 
The influence of Fq shows up through a phase shift, whose wavelength is 0(1): 

A(t, y) = a (y) exp [ -i / F (r, x(r, y)) dr 
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Back to a, we find 

a(t,x)= 7== a(y{t,x))exp(-i [ F (T,x(T,y(t,x))) dr 
\JJt(y(t,x)) V Jo 

where y(t, x) stands for the inverse mapping of y i— > x(t, y). Wigner measures ignore 
this integral, since it corresponds to a phase whose wavelength is large compared 
to e. 

The above approach is very general, and includes linear problems. For instance, 
we may takeF e {x,t) = f(x) £ S(R d ). We could also consider perturbations of order 
0{e a ) with < a ^ 1 instead of 0(e), and follow the same line of reasoning. We 
illustrate this general statement with examples, corresponding to weakly nonlinear 
phenomena. The term "weakly" means that the nonlincarity does not appear in 
the eikonal equation (the geometry of propagation is the same as in the linear case), 
but is present in the leading order transport equation. 

Consider the nonlinear Schrodingcr equation in M. d 

(2.8) ied t r = -^-Ar + V(t,x)r + ef(m 2 )i> s 5 V' e L =0 = aoe^ o/e , 



where / £ C°° (M+ ; R) and Assumption 2.5 is satisfied. Here, 

is a nonlinear function of ip £ . However, we see from (2.7) that, at leading order, 
the modulus of i/j £ is independent of time, so that (2.7) turns out to be a linear 
ordinary differential equation. In [12], the asymptotic behavior of ip 6 is given for 
ie[0,T],by: 

Proposition 2.7 (from [12]). Let f £ C°°(]R + ;R) 7 a £ S(M. d ), and let 
Assumption 2.5 be satisfied. Then for all e £]0, 1], (2.8) has a unique solution 
xjj e £ C°°([0,T] x R d ) n C([0,T];H s ) for all s > d/2 (T is given by Lemma 2.6). 
Moreover, 

V> e - ae iG e iM£ ^0 as e -> 0. 

L oo ([0,T];L 2 nL°°) 

The functions a and G are given by 

X ) = / T 7=7! ' 

G(t, x) = - J f (j,(y(t, ^ |ao(y(*. ^))| 2 ) ds - 

In particular, the unique Wigner measure for ip E is given by: 

w°(t,x,£) = - 1 \a (y(t,x))\ 2 dx<E> ^=v^ elk (t,x)- 

It is independent of the nonlinearity /, and therefore, does not take the nonlinear 
effect, causing the non-trivial presence of G, into account. 

A similar analysis could be carried out by replacing the above local nonlinearity 
by a non-local Hartree type term. We shall simply exhibit an explicit example in 
such a framework. In [15], the following Hartree equation with a harmonic potential 
was considered: 

(2.9) ied t r=-jA^ + ^r+e(\x\-^*m 2 )^ I Vf t=0 = «o, 
with < 7 < 1 and x £ M. d for d ^ 2. The solution to (2.6) is explicit in this case: 

Id 2 7T 

(f> el k(t,x) = — tani, t^-+nZ. 
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Proposition 2.8 ([15, Prop. 4.1]). Let d > 2, a e S(M. d ), and < 7 < 1. 
Let ip £ be the solution to (2.9). Define (for any t) 

-(|x|-^|«o| 2 ) (x) ' 



dT 



II 



I cost|t 

• For ^ t < tt/2, the following asymptotic relation holds: 



sup 



■>P 6 (t,x) 



1 



(cos t)"/ : 



,-«() 



x 



COS T 



•0. 



For tt/2 < t ^ 7T, 



sup 

t<T<7T 



■0 e (r,x) 



(cost)™/ 2 



"0 



(— ) 

V COS T / 



r+jg( 7 



0. 



Por an?/ time t £ [0, 7r]\{f }, the Wigner measure w° associated to the (pure) family 
(ip E )o<e^l is given by: 



1 



cost 



a 



f— ) 

V cos i / 



(ix g) S f= - Xtan t- 



As mentioned in [15], the asymptotic description could be pursued to any time. 
Like in the first example, the Wigner measure is the same as in the linear case and 
ignores leading order nonlinear effects measured by g. 

The phenomenon we described in this paragraph is also present in the main 
result of [14], where a weakly nonlinear perturbation of (1-11) is considered: 

2 

ied t ^ = - E — AiP £ + V(t, xW + V r (-) V> £ + e W\° ?p £ , 

where V is as above, Vr is lattice-periodic, and a £ N. Before the formation 
of caustics, nonlinear effects show up at leading order through a self-modulation 
(phase shift), which may be viewed in this case as a nonlinear Berry phase. 

We emphasize the fact that the self-modulation described in this paragraph 
is not bound to the Schrodingcr equation. In [17], a nonlinear wave equation is 
considered in a weakly nonlinear regime. If in [17], we consider a purely imaginary 
coupling constant, that is, a nonlinear wave equation of the form 

dfu- Au + i\d t u\ p - 1 dtu = 0, 

then we meet the same phenomenon as in this paragraph if, following the notations 
of [17], P_ = or P + = (this corresponds to polarized initial data). 

2.4. Supercritical WKB regime. To conclude on the case of scalar equa- 
tions, consider the case 



(2.10) 



id t ip £ 



m 2 r 



^ £ (0,x) = ag(x)e^°^ )/e . 



We consider the case of a cubic, defocusing nonlincarity for simplicity; the approach 
recalled below can be extended to a wider class a nonlinearities, from [2, 67]. 
Assume that the initial amplitude a§ has an asymptotic expansion of the form 

flo = fl o + ea i + O (e 2 ) , as e — > 0, 

where the functions oq and a\ arc independent of e. The case considered in §2.3 was 
critical as far as WKB analysis is concerned: nonlinear effects are present in the 
transport equation (which determines the leading order amplitude), and it would 
not be the case if the power e in front of the nonlinear term was replaced by e K , 
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with k > 1 (see [12]). The above equation is supercritical as far as WKB analysis 
is concerned. If we seek 

then plugging this expression into (2.10) and ordering the powers of e yields: 
O(e ): 9 t 0+i|V0| 2 + |a o | 2 =O ; 
O (e 1 ) : d t a + V</> • Va + ^aoA</> = 2iRe (a aT) a . 

We see that there is a strong coupling between the phase and the main amplitude: 
ao is present in the equation for <fr. Moreover, the above system is not closed: <fi is 
determined in function of ao, and ao is determined in function of ai. Even if we 
pursued the cascade of equations, this phenomenon would remain: no matter how 
many terms are computed, the system is never closed (see [39]). This is a typical 
feature of supercritical cases in nonlinear geometrical optics (see [20, 21]). 

Suppose however that we know </>, and that the rays associated to <\> do not 
form an envelope for t £ [0, T] (we prefer not to speak of caustic in that case; see 
[10]). Then along these rays, and following the same approach as in 2.3, we see 
that the equation for ao is of the form 

Dt&o = 2iRe (a aT) a . 

Therefore, the modulus of ao is constant along rays, and the coupling between ao 
and ai is present only through a phase modulation for ao. As in 2.3, this phase 
modulation is not trivial in general: a perturbation of the initial data at order 0(e) 
in (2.10) leads to a modification of i\) e at leading order 0(1) for positive times. 
This phenomenon was called ghost effect in a slightly different context [65]. This 
discussion is made rigorous in [12], after rewriting the original idea of E. Grenier 
[45]. Note also that in the above discussion, we have seen that we do not need to 
know ai to determine |ao| 2 : this is strongly related to the following observation. 
Set (p,v) = (|a | 2 , V0). Then the above system for <j>, ao and ai implies 

d t v + v- Vu + V/5 = ; f | t =o = V0 O , 
d t p + V ■(pv) = ; P| t=0 Ha | 2 - 

This system is a polytropic compressible Euler equation. WKB analysis can be 
justified so long as the solution to this system remains smooth ([45, 12]), and the 
Wigner measure is given by 

w°(t, x, £) = p(t, x)dx ® 5 6 = v ( tiX ). 

See also [71]. 

The above remarks can be applied also when ai depends on e. Formally, if 
we replace cl\ by e~ s ai, < S < 1, then we should at least replace ai by e _s ai, 
and also reconsider the remainder 0(e 2 ). Forget this last point. Mimicking the 
above discussion, the interaction between ao and ai leads to a phase modulation 
for ao, of order e~ s : "rapid" oscillations appear. However, such oscillations are 
not detected at the level of Wigner measures: the x-component of the Wigner 
measure ignores this modulation, because it is of modulus one, and the ^-component 
cannot see it, because its wavelength e 5 is too large compared to e. However, this 
phenomenon is everything but negligible at the level of the wave functions, since it 
causes instabilities: 
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Theorem 2.9 (from [11]). Let d^l, a , ai £ S(R d ), (j> G C°°(E rf ;R) ; where 
ao, a± and </>o are independent of e, and V0o G H s (M. d ) for every s ^ 0. Let u £ 
and v £ solve the initial value problems: 

F 2 

Z 3 t u e + — Au £ = \u £ \ 2 u £ ; w e L „ = a e^ o/e . 



it=o 



e 2 



ied t v e + -Av £ = |« e |V ; v £ \ t=Q = (a + e k ai ) e^°' £ , 
for some < k < 1. Assume that Re(aoai) ^ 0. Then we can find < t e — >0 

liminf (t s )-V s (t e )|| L2nioo >0. 
More precisely, this mechanism occurs as soon as t e > e 1 ~ k . Ln particular, 

~ v£ \\l^([0^];L 2 



u \t=o v \t=o 



+oo as e — > 0. 



L 2 

As pointed out above, the Wigner measures ignore this instability. One could 
argue that this instability mechanism does not affect the quadratic quantities, and 
thus may not be physically relevant (the same remark could be made about §2.3). 
Note however that this phenomenon occurs for very small times, when WKB anal- 
ysis is still valid. When WKB ceases to be valid (in this case, this corresponds to 
the appearance of singularities in Euler equations), the approach must be modified. 
We have seen in §2.2 that even for a weaker nonlinearity, nonlinear effects can alter 
drastically the Wigner measures when rays of geometric optics form an envelope. 
The consequence shown in §2.2 was an ill-posedness result for the propagation of 
Wigner measures. Even if Theorem 2.9 may not seem relevant for Wigner mea- 
sures, it may very well happen that for larger times, it causes another ill-posedness 
phenomenon. 

The same discussion remains valid in the case of the Schrodingcr-Poisson sys- 
tem (1.14). It was proven in [72] that so long as the solution to a corresponding 
Euler-Poisson remains smooth, it yields the Wigner measure associated to ip E . 
However, in [1], the approach of Grenier was adapted to the case of (1-14): the 
above instability mechanism can be inferred in this case as well, and the discussion 
remains the same. 

It turns out that this instability mechanism can be met in the case of a weaker 
(as far as WKB analysis is concerned) nonlinearity. Using a semi-classical conformal 
transform, we infer from Theorem 2.9: 

Corollary 2.10 (from [11]). Let d ^ 2, 1 < a < d, and a ,ai G S(R d ) 
independent of e. Let u £ and v e solve the initial value problems: 

2 

ied t u £ + £ —Au e = e a \u £ \ 2 u £ ; u £ (0,x) = a (a;)e- l|;E|2/(2e) , 

2 

ied t v £ + £ —Av £ =s a \v £ \ 2 v £ ; v £ (0,x) = (a (x) + e^^a^x)) e" l|a;|2/(2e) . 
Assume that Re(anai) ^ 0. There exist T £ — > 1~ and < t £ — > such that: 

e-+0 £->0 

(2.11) \\ u£ - v£ \\l~>([o,t?];L2)^° ; lim jn f \\ u£ ~ ve \\L^([o.T-+r^-) > °- 

Consider the case d = 2, and compare with the result of §2.2. The assumption 
a < 2 implies that supercritical nonlinear effects occur near the focusing time t = 1 
(see also [6] for a similar result), but not before since a > 1. This case, along with 
Proposition 2.3, illustrates the above discussion: Wigner measures do not capture 
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the above instability mechanism, but for larger times, this instability may affect 
the Wigner measures. 



3. Limitation of the Wigner measures in the case of systems 

We consider now systems of Schrodingcr equations coupled by a matrix-valued 
potentials. Such systems appear in the framework of Born-Oppcnheimcr approxi- 
mation, where the dynamics of molecules can approximately be reduced to matrix- 
valued Schrodingcr equations on the nucleonic configuration space. We consider 

ied t r = - y AV £ + V(xW, (t, x) g K+ x R d , ^ & C J , 

where the small parameter e is the square root of the ratio of the electronic mass 
on the average mass of molecule's nuclear. Of typical interest is the situation where 



d = 2, j = 2, V(x) 



X\ x 2 
x 2 -x\ 



which is referred to as a codimension 2 crossing (see [47]). Indeed, the eigenvalues 
of V are ±|x|, and they cross on a codimension 2 subspace S of M. d , 

S = {x\ — x 2 = 0}. 

The main features of eigenvalue crossings appear on this example which is simple 
enough so that precise computations can be performed. 

According to the analysis of [42] and taking advantage of the fact that the 
eigenvalues are of multiplicity 1 outside the crossing set S, any Wigner measure w° 
of the family ip E (t) is a 2 by 2 matrix of measures which splits into two parts 

w°(t,x,£) = / + ((,i,On+(i) + /-((,i,on-(3;), 

where H + (x) and II - (x) are the spectral projectors associated respectively with the 
eigenvalues +\x\, —\x\ of V{x) and w°' + , w°'~ are scalar non-negative measures. 
Moreover, outside the crossing set, w 0,± propagate along the classical trajectories 
of ± | a; | according to 

M 

Let us focus on initial data which have only one Wigner measure w® and which 
are microlocally localized on two points of the phase space, so that 

(3.2) o= E 4*(x -4)® 

ie{+,-} 

By [42], as long as the classical trajectories (x ± (t) , ^ (t)) , that is the curves such 
that 

'flki ± (t)=£ ± (t) ; !*(<>) =x±, 

do not reach the crossing set S, ip e {t) has only one Wigner measure w°(t, x,£) which 
is of the form 

w°(t) = a+ (t) 6(x - x+ (t)) ® 6(£ - (*)) n+ + a~ (t) 5{x - aT (t)) <g> 6(£ - C (*)) n ~ 
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with a ± (t) = as long as the trajectories have not reached the crossing. Such 
situations are precisely studied in [34]: it is proved in Proposition 1 of [34] that 
under the assumptions 



4 A e ± = ; ICo" I 2 > 2|aco | ; £o" ■ ^ < 0, 



£o + ' T%T + \/l4 + l 2 + 2141 = -€o" ' 7^7 - V&T " 2|zo I := **, 



the curves (a: (£)) reach S at the same time i* with a non-zero speed £ (f*) 



and one can choose Xq, £q so that they reach S at the same point (0, £*) with 



£* 7^ 0. One then has 

x + (t*) = x-(t*) = 0, t + (t*)=C(t*) 
The trajectories are given on [0, t*] by 



2 \xl\ \xA\ 



2 Fo I Fo I 

Then, there exists ti > t* such that during (t*,ti), both trajectories do not meet 
S again. We will suppose during all this section that we are in this situation and 
for simplicity, we assume \xq\ = 1. We are concerned on the Wigner measure of 
ip £ (t) for t G (t*,ti) which cannot be described by the transport equations above 
mentioned. 

The work [32] proves that the sole knowledge of the initial Wigner measure 
is not enough to determine the Wigner measure of the solution after the crossing 
time t* . A second level of observation is required, and the decisive fact is the way 
the data concentrates on the classical trajectories entering in S with respect to the 
scale y/e. Since all these trajectories are included in the set 

J = {x A £ := xi£ 2 - x 2 £,i = 0}, 

it is enough to study the concentration of ip e {t) on J. 

We come up with the first aspect of our purpose, which is that the sole microlo- 
calization consisting in working in the phase space is not enough. The complete 
analysis of ip E (t) after a crossing point requires a second microlocalization: we add 
to the phase space variables (x, £) a new variable 

i]£l := 1U {+oo, -oo} 

which describes the spread of the wave packet on both sides of J with respect to the 
scale y/e. Then, one defines two-scale Wigner measures whose projections on the 
phase space are the Wigner measure. These measures have been first introduced in 
[60] and developed in [32]; the reader may also find in [30] a survey on the topic. 

In the following, we explain the resolution of our problem by means of two-scale 
Wigner measures (according to [32]). This explains why initial data may have the 
same Wigner measure and generate solutions of the Schrodinger equation which 
have different Wigner measures after a crossing time. We give examples of this 
fact. Then, in a second section, we shall explain and illustrate the limits of the 
two-scale Wigner measures approach: some situations cannot be studied with the 
sole knowledge of the two-scale Wigner measures. 
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3.1. When one needs a second microlocalization. One defines a two- 
scale Wigner measure w°'( 2 ' (t, x, £, 77) of ip £ (t) associated with J, as a weak limit in 
V of the two-scale Wigner functional defined onRj x (M? \ {0}) x R, ; by 



W{t)] (x,t, v ) = w £ W(t)] (x,t)®sL - 

that we test against smooth functions a(x, £,77) which are compactly supported 
in uniformly with respect to 77, and coincide for r\ big enough with an ho- 

mogeneous function of degree in 77, denoted by aoo(a;, £, 77). The knowledge of 
w°'( 2 ' (t, x, £, 77) determines the Wigner measures of ip e (t) above J according to 

Va g C£° (R 2 X x (R| \ {0})) , (w lj,a)=f a(x,£)dw°'W. 

JRJx(R|\{0})xK^ 

Let us suppose that the data ipo nas a unique two-scale Wigner measure w®'^ of 
the form 

w (2) (x,e,7/)= 53 ^-4)®*(f-eg)®4.(»7) ni ^). 

where ^(77) are scalar positive Radon measures on R, and Xq , £^ are as before. 
Define 

(3.3) T{r{) = e-^l"! 2 /!?*! 3 . 

Then, Theorems 2 and 3 in [32] yield 

Proposition 3.1. For any t G [0, f ) U (t*,ti), V> e (£) /ias a unique two-scale 
Wigner measure 77j°'( 2 )(i) which satisfies: 

• Forte [0,t*)> 

7« '^(t,x,c,77)= X! ^-^'W)®^-fW)®^(i)n 3 W- 

• Forte 

ie{+,-} 

Moreover if vf n and v~ n are singular on {\rj\ < 00}, the link between the incident 
measures (v^, v^) and the outgoing ones [ytuti v ~out) * s given by 

1 - T(r?) T(rf) \ (u+ 
T(rj) 1 - T{ V ) J \v- n 

Let us examine the consequence of this Proposition for initial data of the form 
^ = ^ (^P) CX P (i r o + M 2 ) E + (x), I G {1,2,3}, 

with /3i < - = & < Aj, 

where £q~ = r^xj, $ 6 C*5°(R 2 ), and -B + is a smooth bounded function such that 
H + E + = E + on the support of ipn t , and |ji? + ||c2 = 1. These three families have 
the same Wigner measure w® of the form (3.2) with 

4 = \\nh, % =0. 

Simple calculus (see [34, Section 3.2] for details) shows that the two-scale Wigner 
measures associated with J for each of these families are different. Let us denote 
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by uP/^P the two-scale Wigner measure of ip^ £ , £ e {1,2,3}. If v = {vi,v 2 ), we 
denote by v 1 - the vector v 1 = (— v 2 , t>i). Then, we have 

,0,(2), 



w i',i 

0,(2) 



$112 



L- 



5(x-x+)®5(Z-Z+)®5(r)) n+(x), 



W 



1.2 

0,(2) 



/,3 



(x, e, rf) = 6(x - x+) <g> S((, - <g> 1^ R7 (??) ^ n + (x), 
(x,e,r?) =<5(x-x+) <g> (7-^(77-00)+ 7-^(77 + 00)) n+(x), 



with 



' "''O ' 



7(77) 



^(rx+ + r,(x+) ± ) 



dr 



7 



7(77)^77. 



As a consequence of the propagation of two-scale Wigner measures along the clas- 
sical trajectories (see [34]), we get different behaviors after the crossing time t*. 

Corollary 3.2. For £ e {1,2,3}, the family ipf(t) solution to (3.1) with the 
initial data ip^ t has a unique Wigner measure (t) such that 

w°(i, = a+(t)S(x-x + (t))®5(Z-Z + (t))Ti + (x) 

+ aj (t) S(x ~ x- (f)) ® <5(£ - r (*)) IT" (x). 

Moreover, for all t € [0,i*) 

a+(t) = a+(t) = a+(t) = ||$j| 2 L2 , a^(t) = a^(t) = a^(t) = 0, 
and for all t € {t* ,t\), 

4(t) = wnh, or(*) = o, 

a+(t) = (1 - T(7?))7(77), <h(t) = Tfofrfo), 
a+(t)=0, 03W = ||*||| 3 . 

In the first case, the mass propagates along the classical trajectory associated 
with the mode +, while in the third case the mass switches from the mode + to 
the mode — . In the second case, the mass parts between both modes + and — . 

3.2. When the quadratic approach fails. The assumption of singularity 
on the incident measures is crucial in Proposition 3.1. Indeed, if it fails, the two- 
scale Wigner measures after the crossing point cannot be calculated in terms of the 
incident ones. Our aim is now to explain and to illustrate that similar Wigner mea- 
sures (and even similar two-scale Wigner measures) can generate different measures 
after the crossing point. The example given here is precisely discussed in [31]. 

We consider 



r 



X 



X 



,+ 



e- 3 
x — x, 



i Ye 
1 



E UJn 



E+(x) 
E~(x), 



where 



0< 0^0/2 ; 0</3<l/2 



and $ and \& are smooth, compactly supported, functions on M. d , E + (x) (resp. 



E (x)) is a smooth bounded function such that on the support of $ 



(resp. 



of <J 



) one has Il ± (x)E ± (x) = E ± (x) and ||S ± (a;)|| C 2 = 1. We shall focus 



on both situations lUq ^ uj q and ojq = uj q . We set 



= m 



= \m 
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and we suppose 

Vo -=-4(4 A^)^0. 

Proposition 3.3. The family tf> e (t) solution to (3.1) with the initial data 
has a unique Wigner measure w°(t, x,£) such that 

• Forte [0,t*), 

ie{+,-} 

• Forte (t*,h), 

w °(t,x,£)= c>' out 5(x-x j (t)) <8><f IF (i), 

where the coefficients a 3 ' depend on the position of a + and oT with respect to 
1/2: 





+ ,out 


g— ,out 


a-,a + < 1/2 
a+ < oT = 1/2 
a~ < a + = 1/2 
f a+ = oT = 1/2 

1 Vo ^ Vo 


c +,in 

c +,in +T(r)o)c-> in 

c +,in (1 _ T{r] +^ 
c +,in (1 _ r(?? +)) 

+T(r /0 -)c-' i " 


c -,tn 

C --»(1-T( % -)) 

c -,in + T(r/+ )c+' m 

c -^(l-T( v -)) 
+T(r / +)c+' i ™ 



If a + = a =1/2 and ?7q = r/ := 970, there exists po e M + and 4>o e K smc/j i/iai 
c+>out = c+ , in (1 _ T(??o)) + r(r?o)c -,™ + pQ cos(0o); 

c -,o« t = c -,™ (1 _ T(r?())) + T(r?0 ) c +,™ _ pQ cos (^ ). 

Besides, for any <j> € M, if one turns ^ into e^'I' , t/ien one /las 

c+>0 «t = c+!i „ (1 _ T(??o)) + T(??0 ) c -,m + po cos(0o _ ^ 

c _, out = c _, m (1 _ T(??o)) + T{m)c+ ,n _ po cos(0Q _ ^ 

Here again, we see that the interaction between both incident modes cannot 
be described only by the knowledge of the Wigner measure: one needs to know the 
two-scale Wigner measure. Indeed, for each case described in the above array, the 
two-scale Wigner measure is different: 

ie{+,-} 

with v 3 (rf) supported on 00 for a 3 < 1/2 and v 3 (vf) supported on ?/q if a 3 = 1/2. 
Therefore, the situations of the array enters in the range of validity of Proposi- 
tion 3.1. In the last situation where a + = a~ =1/2 and 7]q = 7]q , the two incident 
wave functions interact and some coupling term appears. By modifying the initial 
data, turning ^ into e lX ^, one does not modify the two-scale Wigner measure but 
one can transform the coupling (and even suppress it) so that the outgoing Wigner 
measures change. The same thing happens for any initial data provided it has the 
same two-scale Wigner measure as our example. 

For proving Proposition 3.3, we use the normal form one can obtain for (3.1), 
and which is crucial for proving Proposition 3.1. The reader can refer to [23] for 
the more elaborate result for general systems presenting codimension 2 crossings 
(see also [24] for the case of hermitian matrix- valued symbol). One can find in [32] 
a weaker result which covers the case of Schrodinger equation (3.1), and is enough 
for calculating Wigner measures. Through a change of symplectic coordinates in 
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space-time phase space, and a change of unknown by use of a Fourier Integral 
operator, one reduces to the system: 



(3.5) 



: d s U E 



Zl 



(s,z = {z 1 ,z 2 )) e 



We outline two important features of this normal form. On the one hand, this 
normal form is microlocal, near the crossing point (t* , 0, t*. £*) where r* is the 
If *l 2 

energy variable r* = . On the other hand, space time coordinates are involved 
by the canonical transform (i.e. by the change of symplcctic coordinates) 

K : (t,X,T,£) i * (s,z,a,(). 

Through n, the crossing set S becomes S = {s = Z\ = 0} and the set J = {x A £ = 
0} becomes J = {z\ = 0}. Moreover, the classical trajectories in space time 
coordinates (r,x ± (r),T,£ ± (r)) r>0 (where r = — i|£ ± (r)| 2 ±|x ± (r)| = Const.) which 
enters in S maps on the curves (r, 0, z%, ±|r|, Ci, C2) r >o an d the link between Z\ and 
x A £ is given by Z\ = e(t, x, r, £) x A £, where e is a smooth function such that 
e(t*, 0, r*, £*) = |£*| -3 / 2 . Finally, the Wigner measure is invariant by this change 
of coordinates, and if we denote by A the additional variable of the two-scale Wigner 
measure in the variables (s, z), and by r\ the corresponding variables in (x, £), we 
have A = e(t, x, r, because of the link between z\ and (t, £, r, £). 

We focus now on system (3.5). This system of o.d.e. is simple enough so that 
direct calculations are possible. This has been done by Landau and Zener in the 
30's (see [52] and [70]). We use the description of u e near Z\ = as stated in [32, 
Proposition 9], which is obtained by means of stationary phase method. One can 
also find a resolution of this system in [47], where special functions are used. 

PROPOSITION 3.4. There exist families of vectors ofC 2 , a e = (a\,a\), u E = 
(u^wf), such that, as e goes to and for bounded, 
• For s < 0, 



u\(s,z) 
For s > 0, 



«1 



n 9 



o(l). 



u\(s,z) 



Moreover 



o(l), 



(3.6) S(X) 



o(A) -6(A) 
6(A) <z(A) 



with 



a(A) 



o(l) 



a(A) 2 + |6(A)| 2 = l. 



Proof of Proposition 3.3. Let us use the relation between the classical tra- 
jectories to identify v + ' m (resp. is~' m ) as the two-scale Wigner measure of a\ 
(resp. af) f° r { z i = 0}, and i/+> ou * (resp j/ - > out ) as the one of wf (resp. wf). 
• If (a* < = 1/2) or (a* = 1/2 and 7^J~ 7^ "Ho), the two incoming measures are 
mutually singular. Therefore, the measure of 

u\(z) = a f-^Lj a\(z) - b f-^=J a £ 2 (z) 

is the sum of the measures of each term. We obtain 

v-° ut (z, C, A) = a(A) 2 ^"(z, C, A) + \b{\)\ 2 v- m {z , C, A), 
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which gives in the variables (x, £, 77) 

v-'° ut (x, t, V ) = T{rf)v + ' in {x, & 77) + (1 - Tfa)) v~' in (x, £, V ), 
where we have used 

a (v irr 3/2 ) 2 = t{ V ), \b (r? irr 3/2 ) f = i - m. 

Similarly, we get 

v + '° ut (x, £, r,) = (1 - T( V )) v + - m {x, t rj) + T{r,)u-' in (x, f , ?7 ). 

We see here how the singularity relation plays a role for finite 7/ in Proposition 3.1. 

• If or* 1 < 1/2, both incident measures are localized in r\ = oo, Prop. 3.1 yields 

v ± >° ut {x,i,ri)=v ± > in {x,li,n)- 

• If = 1/2 and t]q = t]q . Observing that u ± ' out is localized above = t]q := T)o 
and taking into account the two-scale joint measure 9 between a\ and a\ which is 
also supported above 770 , we obtain 

= (1 - T(r, )) v+> m + T(r, )p-' in + 2 Re (a \C r 3/2 ) b (r? \C T 3/2 ) fl) , 

T(7y )^+^ m + (1 - T(r /0 )) v~> m - 2 Re (a (770 |ri" 3/2 ) 6 (jla \C\ 3/2 ) 0) ■ 

There exists (po,<po) £ K + x [0, 2tt] such that 

«(^oirr 3/2 )&(%irr 3/2 ) = P0 e i *>. 

One then observes that if one multiplies the minus component of the initial data 
by e 4 "^, one turns a| into e^a| an d ^ into e~ l< ^9, whence the result. □ 



v 



[1 

[2, 
[3 
['1 
[5 
[6 

[7: 

[8 
[9 
[10 

[11 
[12 
[13 

[11 
[15 



References 

T. Alazard and R. Carles, Semi- classical limit of Schrddinger-Poisson equations in space 
dimension n ^ 3, J. Differential Equations 233 (2007), no. 1, 241-275. 

T. Alazard and R. Carles, Supercritical geometric optics for nonlinear Schrddinger equations, 
archived as arXiv: 0704. 2488, 2007. 

P. Bechouche, N. J. Mauser, and F. Poupaud, Semiclassical limit for the Schrddinger equation 
in a crystal with Coulomb interaction, Comm. Pure Appl. Math. 54 (2001), no. 4, 851-890. 
N. Burq, Mesures semi-classiques et mesures de defaut, Asterisque (1997), no. 245, Exp. No. 
826, 4, 167-195, Seminaire Bourbaki, Vol. 1996/97. 

N. Burq, Semiclassical estimates for the resolvent in non trapping geometry, Int. Math. Res. 
Notices 5 (2002), 221-241. 

N. Burq and M. Zworski, Instability for the semiclassical non-linear Schrddinger equation, 
Comm. Math. Phys. 260 (2005), no. 1, 45-58. 

R. Carles, Geometric optics with caustic crossing for some nonlinear Schrddinger equations, 
Indiana Univ. Math. J. 49 (2000), no. 2, 475-551. 

R. Carles, Remarques sur les mesures de Wigner, C. R. Acad. Sci. Paris, t. 332, Seric I 332 
(2001), no. 11, 981-984. 

R. Carles, Semi-classical Schrddinger equations with harmonic potential and nonlinear per- 
turbation, Ann. Inst. H. Poincare Anal. Non Lineairc 20 (2003), no. 3, 501-542. 
R. Carles, Cascade of phase shifts for nonlinear Schrddinger equations, J. Hyperbolic Differ. 
Equ. 4 (2007), no. 2, 207-231. 

R. Carles, Geometric optics and instability for semi- classical Schrddinger equations, Arch. 
Ration. Mech. Anal. 183 (2007), no. 3, 525-553. 

R. Carles, WKB analysis for nonlinear Schrddinger equations with potential, Comm. Math. 
Phys. 269 (2007), no. 1, 195-221. 

R. Carles, Semi-classical analysis for nonlinear Schrddinger equations, World Scientific Pub- 
lishing, 2008. 

R. Carles, P. A. Markowich, and C. Sparber, Semiclassical asymptotics for weakly nonlinear 
Bloch waves, J. Stat. Phys. 117 (2004), no. 1-2, 343-375. 

R. Carles, N. J. Mauser, and H. P. Stimming, (Semi)classical limit of the Hartree equation 
with harmonic potential, SIAM J. Appl. Math. 66 (2005), no. 1, 29-56. 



24 R. CARLES, C. FERMANIAN, N. MAUSER, AND H. P. SUMMING 

[16] R. Carles and T. Ozawa, On the wave operators for the critical nonlinear Schrodinger equa- 
tion critical nonlinear Schrodinger equation, Math. Res. Lett. 15 (2008), no. 1, 185-195. 

[17] R. Carles and J. Rauch, Focusing of spherical nonlinear pulses in R 1 + 3 , Proc. Amer. Math. 
Soc. 130 (2002), no. 3, 791-804. 

[18] F. Castella, Propagation of space moments in the Vlasov-Poisson equation and further re- 
sults, Ann. Inst. H. Poincare Anal. Non Lineaire 16 (1999), no. 4, 503—533. 

[19] T. Cazenave, Semilinear Schrodinger equations, Courant Lecture Notes in Mathematics, 
vol. 10, New York University Courant Institute of Mathematical Sciences, New York, 2003. 

[20] C. Cheverry, Cascade of phases in turbulent flows, Bull. Soc. Math. France 134 (2006), no. 1, 
33-82. 

[21] C. Cheverry and O. Cues, Counter-examples to concentration- cancellation and supercritical 
nonlinear geometric optics for the incompressible Euler equations, Arch. Ration. Mech. Anal. 
(2007), to appear. 

[22] Y. Colin de Verdiere, Ergodicite et fonctions propres du laplacien, Comm. Math. Phys. 102 
(1985), no. 3, 497-502. 

[23] Y. Colin de Verdiere, The level crossing problem in semi- classical analysis. I. The symmetric 

case, Ann. Inst. Fourier (Grenoble) 53 (2003), no. 4, 1023-1054. 
[24] Y. Colin de Verdiere, The level crossing problem in semi-classical analysis. II. The Hermitian 

case, Ann. Inst. Fourier (Grenoble) 54 (2004), no. 5, 1423-1441, xv, xx-xxi. 
[25] J. Derezihski and C. Gerard, Scattering theory of quantum and classical N-particle systems, 

Texts and Monographs in Physics, Springer Verlag, Berlin Heidelberg, 1997. 
[26] R. DiPerna and P.-L. Lions, Solutions globales d'equations du type Vlasov-Poisson, C. R. 

Acad. Sci. Paris Scr. I Math. 307 (1988), no. 12, 655-658. 
[27] W. Domcke, D. Yarkony, and H. Koppel, Conical intersections, Worl Scientific Publishing, 

2004. 

[28] C. Fermanian Kammerer, Propagation and absorption of concentration effects near shock 

hypersurfaces for the heat equation, Asymptot. Anal. 24 (2000), no. 2, 107-141. 
[29] C. Fermanian Kammerer, Semiclassical analysis of generic codimension 3 crossings, Int. 

Math. Res. Not. 45 (2004), 2391-2435. 
[30] C. Fermanian Kammerer, Analyse a deux echelles d'une suite bornee de L 2 sur une sous- 

variete du cotangent., C. R. Acad. Sci. Paris, t. 332, Seric I 340 (2005), 269-274. 
[31] C. Fermanian Kammerer, Normal forms for conical intersections in quantum chemistry, 

Mathematical Physics Electronic Journal (2007). 
[32] C. Fermanian Kammerer and P. Gerard, Mesures semi-classiques et croisements de modes, 

Bull. Soc. math. France 130 (2002), no. 1, 123-168. 
[33] C. Fermanian Kammerer and P. Gerard, A landau-zener formula for non- degenerated invo- 

lutive codimension 3 crossings, Ann. Henri Poincare 4 (2003), 513-552. 
[34] C. Fermanian Kammerer and C. Lasser, Wigner measures and codimension two crossings, 

Jour. Math. Phys. 44 (2003), no. 2, 507-527. 
[35] C. Fermanian Kammerer and V. Rousse, Resolvant estimates and matrix-valued Schrodinger 

operator with eigenvalue crossings. Application to Strichartz inequalities, Comm. in Partial 

Dili. Eq. 33 (2008), no. 1, 19-44. 
[36] I. Gasser, P.-E. Jabin, and B. Perthame, Regularity and propagation of moments in some 

nonlinear Vlasov systems, Proc. Roy. Soc. Edinburgh Sect. A 130 (2000), no. 6, 1259-1273. 
[37] P. Gerard, Mesures semi-classique et ondes de Bloch, Seminaire sur les Equations aux 

Derivees Partielles, 1990-1991 (Palaiseau), Ecolc Polytech., 1991, pp. Exp. No. XVI, 19. 
[38] P. Gerard, Microlocal defect measures, Comm. Part. Diff. Eq. 16 (1991), 1761-1794. 
[39] P. Gerard, Remarques sur I'analyse semi-classique de I'equation de Schrodinger non lineaire, 

Seminaire sur les Equations aux Derivees Partielles, 1992-1993, Ecole Polytech., Palaiseau, 

1993, pp. Exp. No. XIII, 13. 
[40] P. Gerard, Oscillations and concentration effects in semilinear dispersive wave equations, J. 

Funct. Anal. 141 (1996), no. 1, 60-98. 
[41] P. Gerard and E. Leichtnam, Ergodic properties of eigenf unctions for the Dirichlet problem, 

Duke Math. J. 71 (1993), no. 2, 559-607. 
[42] P. Gerard, P. A. Markowich, N. J. Mauser, and F. Poupaud, Homogenization limits and 

Wigner transforms, Comm. Pure Appl. Math. 50 (1997), no. 4, 323-379. 
[43] P. Gerard, P. A. Markowich, N. J. Mauser, and F. Poupaud, Erratum : Homogenization 

limits and Wigner transforms, Comm. Pure Appl. Math. 53 (2000), 280-281. 
[44] J. Ginibre and G. Velo, On a class of nonlinear Schrodinger equations with nonlocal inter- 
action, Math. Z. 170 (1980), no. 2, 109-136. 
[45] E. Grenier, Semiclassical limit of the nonlinear Schrodinger equation in small time, Proc. 

Amer. Math. Soc. 126 (1998), no. 2, 523-530. 



WIGNER MEASURES FOR SCHRODINGER EQUATIONS 



25 



[46] A. Grigis and J. Sjostrand, Microlocal analysis for differential operators, London Mathemat- 
ical Society Lecture Note Series, vol. 196, Cambridge University Press, Cambridge, 1994, An 
introduction. 

[47] George A. Hagedorn, Molecular propagation through electron energy level crossings, Mem. 

Amcr. Math. Soc. Ill (1994), no. 536, vi+130. 
[48] N. Hayashi and Y. Tsutsumi, Scattering theory for Hartree type equations, Ann. Inst. H. 

Poincare Phys. Theor. 46 (1987), no. 2, 187-213. 
[49] B. Hclffer, A. Martinez, and D. Robert, Ergodicite et liraite semi-classique, Comm. Math. 

Phys. 109 (1987), no. 2, 313-326. 
[50] T. Jecko, From classical to semiclassical non-trapping behaviour, C. R. Acad. Sci. Paris Ser. 

I Math. 338 (2004), 545-548. 
[51] T. Jecko, Non-trapping condition for semiclassical Schrddinger operators with matrix-valued 

potentials, Math. Phys. Electron. J. 11 (2005), Paper 2, 38 pp. 
[52] L. Landau, Collected papers of L. Landau, Pergamon Press, 1965. 

[53] G. Lebeau, Equation des ondes amorties, Algebraic and geometric methods in mathematical 
physics (Kacivcli, 1993), Math. Phys. Stud., vol. 19, Kluwer Acad. Publ., Dordrecht, 1996, 
pp. 73-109. 

[54] P.-L. Lions and T. Paul, Sur les mesures de Wigner, Rev. Mat. Iberoamericana 9 (1993), 
no. 3, 553-618. 

[55] P.-L. Lions and B. Perthame, Propagation of moments and regularity for the 3- dimensional 
Vlasov-Poisson system, Invent. Math. 105 (1991), no. 2, 415-430. 

[56] A. J. Majda, G. Majda, and Y. X. Zheng, Concentrations in the one- dimensional Vlasov- 
Poisson equations. I. Temporal development and non-unique weak solutions in the single 
component case, Phys. D 74 (1994), no. 3-4, 268-300. 

[57] P. A. Markowich and N. J. Mauser, The classical limit of a self- consistent quantum- Vlasov 
equation in 3D, Math. Models Methods Appl. Sci. 3 (1993), no. 1, 109-124. 

[58] P. A. Markowich, N. J. Mauser, and F. Poupaud, A Wigner-function approach to 
(semi) classical limits: electrons in a periodic potential, J. Math. Phys. 35 (1994), no. 3, 
1066-1094. 

[59] N. J. Mauser, (Semi) classical limits of Schrddinger-Poisson systems via Wigner transforms, 

Journecs "Equations aux Derivees Partielles" (Forges-les-Eaux, 2002), Univ. Nantes, Nantes, 

2002, pp. Exp. No. XI, 12. 
[60] L. Miller, Propagation d'ondes semi-classiques a travers une interface et mesures 2- 

microlocales., Ph.D. thesis, Ecole polytcchnique, 1996. 
[61] F. Nier, Une description semi-classique de la diffusion quantique, Scminairc sur les Equations 

aux Derivees Partielles, 1994-1995, Ecole Polytech., Palaiseau, 1995, pp. Exp. No. VIII, 10. 
[62] F. Nier, A semi-classical picture of quantum scattering, Ann. Sci. Ecole Norm. Sup. (4) 29 

(1996), no. 2, 149-183. 

[63] G. Papanicolau and J. B. Keller, Stochastic differential equations with applications to random 
harmonic oscillators and wave propagation in random media, SIAM J. Appl. Math. 21 (1971), 
287-305. 

[64] A. I. Snircl'man, Ergodic properties of eigenf unctions, Uspehi Mat. Nauk 29 (1974), 
no. 6(180), 181-182. 

[65] Y. Sone, K. Aoki, S. Takata, H. Sugimoto, and A. V. Bobylev, Inappropriateness of the 
heat- conduction equation for description of a temperature field of a stationary gas in the 
continuum limit: examination by asymptotic analysis and numerical computation of the 
Boltzmann equation, Phys. Fluids 8 (1996), no. 2, 628-638. 

[66] L. Tartar, H-measures, a new approach for studying homogenisation, oscillations and con- 
centration effects in partial differential equations, Proc. Roy. Soc. Edinburgh Sect. A 115 
(1990), no. 3-4, 193-230. 

[67] L. Thomann, Instabilities for supercritical Schrddinger equations in analytic manifolds, J. 
Differential Equations (2008), to appear. Archived as arXiv : 0707 . 1785. 

[68] E. Wigner, On the quantum correction for thermodynamic equilibrium, Phys. Rev. (2) 40 
(1932), 749-759. 

[69] S. Zelditch, Uniform distribution of eigenfunctions on compact hyperbolic surfaces, Duke 

Math. J. 55 (1987), no. 4, 919-941. 
[70] C. Zener, Non-adiabatic crossing of energy levels, Proc. Roy. Soc. Lond. 137 (1932), 696-702. 
[71] P. Zhang, Semiclassical limit of nonlinear Schrddinger equation. II, J. Partial Differential 

Equations 15 (2002), no. 2, 83-96. 
[72] P. Zhang, Wigner measure and the semiclassical limit of Schrddinger-Poisson equations, 

SIAM J. Math. Anal. 34 (2002), no. 3, 700-718. 



2<> 



R. CARLES, C. FERMANIAN, N. MAUSER, AND H. P. STIMMING 



[73] P. Zhang, Y. Zheng, and N. J. Mauser, The limit from the Schrddinger-Poisson to the Vlasov- 
Poisson equations with general data in one dimension, Comm. Pure Appl. Math. 55 (2002), 
no. 5, 582-632. 

[74] Y. X. Zheng and A. Majda, Existence of global weak solutions to one-component Vlasov- 
Poisson and Fokker- Planck- Poisson systems in one space dimension with measures as initial 
data, Comm. Pure Appl. Math. 47 (1994), no. 10, 1365-1401. 

[75] P. Zhidkov, On global solutions for the Vlasov-Poisson system, Electron. J. Differential Equa- 
tions (2004), No. 58, 11 pp. (electronic). 

(R. Carles) CNRS & Universite Montpellier 2, Mathematiques, CC 051, Place Eugene 
Bataillon, 34095 Montpellier cedex 5, France 
E-mail address: Remi.CarlesOmath.cnrs.fr 

(C. Fermanian) LAMA UMR CNRS 8050, Universite Paris EST, 61, avenue du General 
de Gaulle, 94010 Creteil Cedex, France 

E-mail address: Clotilde.Fermanianauniv-parisl2.fr 

(N.J. Mauser) Wolfgang Pauli Institute c/o Fak. f. Math., Univ. Wien, Nord- 
bergstr. 15, A 1090 Wien, Austria 

E-mail address: mauser3courant.nyu.edu 

(H. P. Stimming) Wolfgang Pauli Institute c/o Fak. f. Math., Univ. Wien 
E-mail address: nans .peter . stimmingOunivie . ac . at 



